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ABSTRACT 

Aims. An extreme Kerr black hole (BH) surrounded by a precessing disk is invoked to explain the light curves of 
gamma-ray bursts (GRBs) based on the coexistence of the Blandford-Znajek (BZ) and the magnetic coupling (MC) 
processes. 

Methods. The overall shape of the light curves and the duration of GRBs are interpreted by the evolution of the half- 
opening angle of the magnetic flux on the BH horizon, and the complex temporal structures are modulated by the 
precession and nutation of the jet powered by the BZ process. 

Results. The time profile of the emission exhibits a fast rise and a slow decay due to the effect of the evolution of the 
half-opening angle. The light curves of several GRBs are well fitted by this model with only six free parameters. 

Key words, gamma-ray: bursts - black hole - accretion disk - Blandford-Znajek process 



1. Introduction 

Gamma ray bursts (GRBs) are possibly the most luminous 
objects in the Universe. Extremely high energy released in 
very short timescales suggests that GRBs involve the for- 
mation of a black hole (BH) via a catastrophic stellar col- 
lapse event or possibly a neutron star merger, implying that 
an inner engine could be built on an accreting BH (Piran 
2004). 

Among a variety of mechanisms for powering GRBs, 
the BZ process (Blandford & Znajek 1977) has the unique 
advantage of providing "clean" (free of baryonic contam- 
ination) energy by extracting rotating energy from a BH 
and transferring it in the form of Poynting flow in the out- 
going energy flux (Lee et al. 2000). Recently, observations 
and theoretical considerations have linked long-duration 
GRBs with ultrabright Type Ib/c SNe (Galama et al. 1998; 
Bloom et al. 1999). Brown et al. (2000) proposed a specific 
scenario for a GRB-SN connection. They argued that the 
GRB is powered by the BZ process and the SN is pow- 
ered by the MC process, which is regarded as one of the 
variants of the BZ process (van Putten 1999; Blandford 
1999; Li 2000, 2002; Wang et al. 2002). However, they 
failed to distinguish the fractions of the energy for these 
two objects. More recently, van Putten and collaborators 
proposed a dominant spin-connection between the central 
BH and surrounding high-density matter (van Putten 2001; 
van Putten & Levinson 2003). It is based on similar shapes 
in topology of the torus magnetosphere with the magneto- 
sphere of pulsars, when viewed in a poloidal cross-section. 
This description points towards complete calorimetry on 
GRB-SNe, upon including an unexpectedly powerful long- 
duration burst of gravitational-radiation, representing most 
of the spin-energy liberated from the central BH. 

Lei et al. (2005, hereafter LWM05) proposed a scenario 
for GRBs in Type Ib/c SNe, invoking the coexistence of 
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the BZ and MC processes. In LWM05 the half-opening an- 
gle of the magnetic flux tube on the horizon is determined 
by the mapping relation between the angular coordinate 
on the BH horizon and the radial coordinate on the sur- 
rounding accretion disk. In this scenario the half-opening 
angle evolves to zero with the spinning-down BH. This ef- 
fect shuts off a GRB, and the overall timescale of the GRB 
could be fitted by the lifetime of the open magnetic flux on 
the horizon. 

Besides the feature of high energy released in very short 
durations, most GRBs are characterized by highly vari- 
able light curves with complex temporal behavior. The 
usual explanations of the temporal structures range from 
multiple shock fronts running into the ambient medium 
(Sari, Narayan & Piran 1996), expanding shells with 
brighter patches and dimmer regions (Fenimore, Madras & 
Nayakshin 1996), to repeated series of pulses with Gaussian 
or power-law profiles (Norries et al. 1996). However, a 
clear physical explanation is lacking. Some authors (Fargion 
& Sails 1996; Blackman, Yi & Field 1996; Fargion 1999; 
Portegies Zwart, Lee & Lee 1999, hereafter PZLL99; 
Portegies Zwart & Totani 2001, hereafter PZTOl; Fargion 
& Grossi 2006) suggested that the light curves of GRBs 
can be explained by a beamed emission from a precessing 
jet. PZLL99 constructed ad hoc a function I{t) to describe 
the outline of the light curves, which is characterized by 
three timescales: a fast rise with timescale Trise^ a plateau 
phase with timescale Tpiat , and a stiff decay with timescale 
Tdecay The complex temporal structures are modulated by 
the precession and nutation of the jet. Although their model 
successfully fits several light curves of GRBs, the origin of 
the function I{t) is lacking. 

In fact, the evolution of the half-opening angle on the 
BH horizon described in LWM05 provides a natural expla- 
nation for the outline of the light curves. Based on LWM05 
and PZLL99, in this paper, we intend to combine the evo- 
lution of the half-opening angle with the precession and 
nutation of the jet to fit the complex light curves of GRBs. 
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Fig. 1. The polodial topology of the magnetosphere of a 
disk surrounding a rapidly rotating BH (not to scale). Along 
the open magnetic flux tube, an asymptotic coUimated jet 
may be formed. 

The paper is organized as follows. In section 2 we de- 
rive the jet luminosity per steradian by the mapping rela- 
tion between angular coordinate on the BH horizon and the 
cylindrical radius of the coUimated jet. Based on the evo- 
lution of an extreme Kerr BH we obtain an intrinsic time 
variation of the gamma-ray emission. It is shown that the 
time profile of the emission exhibits a fast rise and a slow 
decay. In section 3, the light curves of several GRBs are 
fitted by combining the evolution of the half-opening angle 
with the precession and nutation of the jet. In section 4, we 
summarize the main results and discuss some issues related 
to our model. Throughout this paper the geometry units 
G = c ~1 are used. 



2. Time profile of jet luminosity 

The poloidal configuration of the magnetic field is shown 
in Fig. 1, which is adapted from van Putten (2001). 

In Fig. 1 the angle Oh is the half-opening angle of the 
open magnetic flux tube, indicating the angular boundary 
between open and closed field lines on the horizon. The 
angle 9h can be determined by (Wang et al., 2003) 



Fig. 2. The curves of Oh vs a, for n=3.5, 4, 4.5, 5, 5.493. 
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Eq. 11} is derived based on the conservation of magnetic 
flux with three assumptions: (i) the magnetic field on the 
horizon Bh is constant, (ii) the magnetic field on the disk 
surface varies as a power law with the radial coordinate 
of the disk, and (iii) the magnetic flux connecting the BH 
with its surrounding disk takes precedence over that con- 
necting the BH to the remote load. In Eq. ([T]) a* = J/AP 
is the BH spin defined in terms of the BH mass M and the 
angular momentum J, the parameter n is the power-law in- 
dex for the variation of Bd, i.e., Bg^"", and ^ = r/vms is 
the radial coordinate on the disk, which is defined in terms 
of the radius Vms = -^Xms "^f the marginally stable orbit 
(Novikov & Thorne 1973). The hmction G{a^;^, n) is given 

by 



G(a,;^,n) 
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By using Eq. ^ we have the curves of 6h versus a* 
with different values of the power-law index n as shown in 
Fig. 2, in which the half-opening angle Oh always increases 
monotonically with increasing n for a given BH spin a*. 

A very interesting feature shown in Fig. 2 is that the an- 
gle 9h evolves non-monotonically with decreasing a*, i.e., it 
increases very rapidly as a* spins down from unity and then 
decreases slowly to zero for 3.003 < n < 5.493. Therefore, 
the time-scale of the duration of the burst is given by the 
lifetime of rapid spin of the BH as mentioned by van Putten 
(2001). This evolution characteristic can be applied to fit- 
ting GRBs in two aspects: ([IJ determining the duration 
of the GRBs as described in LWM05, and ^ shaping the 
time profile of jet luminosity. To see the latter effect, we 
will derive the jet luminosity in §2.1, and investigate its 
time profile in §2.2. 

2.1. Jet luminosity 

As is well known, the rotating energy is extracted from 
a BH in the form of Poynting flow in the BZ process, 
and this energy dissipates and accelerates electrons to pro- 
duce gamma-rays through synchrotron radiation or inverse 
Compton scattering. Although some works approach the 
magnetic dissipation in jets, its origin remains unclear 
(Lyutikov & Blandford 2003; Spruit, Daigne, & Drenkhahn 
2001). As a simple analysis, we introduce a parameter to 
denote the fraction of BZ energy converted into gamma-ray 
energy. 

The BZ power transferred through two adjacent mag- 
netic surfaces between 6 and 6 + d9 on the BH horizon is 
given by (Wang et al. 2002) 



dpBz = dPez/Po = 2al 



fc(l - k) sin^ e 
2-{l-q) sin^e' 



d0. 



(3) 
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where Pq = Bj^AI^, q = -^/l — , and k is the ratio of the 
angular velocity of open magnetic field lines to that of the 
BH. For the optimal BZ power, k = 0.5 is taken. 

There are a variety of approaches to studying coUimat- 
ing jet flows, such as the curvature force exerted by the 
toroidal field (Sakurai 1985), the pressure of the poloidal 
magnetic field surrounding the jet (Blandford 1993; Spruit 
1994), and the pressure of the baryon-rich wind supported 
by the surrounding disk (van Putten 2001; van Putten & 
Levinson 2003). The coUimation of the flow is determined 
by the cross-field force-balance in the direction perpendic- 
ular to the field. Fcndt (1997, hereafter F97) obtained the 
numerical solution of the stream equation for this force- 
balance of a force-free magnetic jet. The asymptotic jet 
can be coUimated into a cylindrical shape with a jet radius 
of several light cylinder radii. 

In order to work out an analytical model we adopt 
F97's scenario and assume that the poloidal component 
of the magnetic field of the coUimated jet, Bl, varies as 
Bl/Bh = e~'^^, where x = R/th is the cylindrical radius 
of the coUimated jet in terms of the radius of the BH hori- 
zon. Based on the conservation of magnetic flux, we have 

= Bjji: {wp)^^^^ dO = BlnRdR. (4) 

Substituting {'cop)^^^^ — X]r=rjf ^i'^^ ~ 2MrH svaO into 
Eq. (HI), we have 

2suiede/{l + q) = xe-^^dx. (5) 

Integrating Eq. (O, we obtain the mapping relation be- 
tween the angular coordinate 6 on the BH horizon and the 
cylindrical radius x the coUimated jet as follows 

cos = 1 - [1 - (1 + Ax)e-^"] . (6) 

Defining Rl = c/flp — 2/51 as the light cylindrical 
radius of the coUimated jet, we have xl = Rl/th — 4/a*. 
Following F97, the radius Rjet of a cylindrically coUimated 
jet should be less than AR^. In this work, we assume Rj^t ~ 
Rl, and the field lines from Oh near the BH will connect 
to Rl in the coUimated region as shown in Fig. 1. Thus the 
relation between 9h and is expressed by 

cos^ff = l-^[l-(l + Aa;L)e-^^'-]. (7) 

Following PZLL99, the angle dependence of the jet is 
given approximately as x = Xfsm^, where xj = rf/rn 
is the freezeout distance from the BH to the remote load, 
where the photons can leave freely. The opening angle ipjet 
is fixed by smipjet ~ Xhjxf as PZLL99, and a small open- 
ing angle ■i/'jet = 6° with Xf = lOa;^ is assumed. 

Based on the above equations we express the jet lumi- 
nosity per steradian for an observer at angle with respect 
to the central locus of the jet, i.e., 

Ljetia^^OH^i!) = Ljet/Po = BydPsz/dfl 

_ EjdPBZ 

27r sin TpdTp 

_ £7 dpBz d9 dx 
27r sin ^ dO dx dtp ' 
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Fig. 3. V-dependence of Ljet a) for a* — 0.99 with n —5, 
4.5, 4 and 3.5 in solid, dash-dotted, dashed and dotted lines, 
respectively, and b) for n = 5 with a* = 0.9999, 0.9, 0.8 and 
0.7 in solid, dash-dotted, dashed and dotted fines, respec- 
tively. 

where 6h = OHia*,n) is determined by Eq. ([T]). 

Taking ~ 0.15 in calculations as given by van Putten 
et al. (2004), we have Ljet varying with tp as shown in Fig. 
3. From Fig. 3 we find that the value of Ljet is sensitive to n 
and a*. It vanishes at ■(/' = 0, and attains its peak value near 
the spin axis of the BH. In addition, we find that both the 
value and the angular coordinate of the peak Ljet strongly 
depend on the values of a* and n, while the variation of 
Ljet with tp remains similar. 

The 6h dependence of Ljet for the given values of a* 
and -ip is depicted in Fig. 4, and we find that the value of 
Ljet is sensitive to that of 6h- 

Combining the evolution of 6h with the BH spin, we 
have the curves of the luminosity Ljet versus a* as shown 
in Fig. 5. 

Inspecting Fig. 5, we find that the luminosity Ljet 
evolves non-monotonically with the spin of an extreme 
Kerr BH, increasing very fast at first and then decreas- 
ing slowly. This result is understandable by considering the 
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Fig. 4. The curves of Ljet vs Oh for a* = 0.9 and ip — 
l°(solid line), a* = 0.7 and ?/> = 1° (dash-dotted line), and 
a* = 0.9 and tp = 3°(dashed line). 



0.5 



0.4 



0.3 



0.2 



0.1 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




0.8 0.85 0.9 0.95 1 



Fig. 5. The curves of Ljet vs a* a) for ip = 1° with n =5, 
4.5, 4 and 3.5 in solid, dash-dotted, dashed and dotted lines, 
respectively; b) for n = 4 with ip = 1°,2°,3° and 4° in solid, 
dash-dotted, dashed and dotted lines, respectively. 



non-monotonic evolution of 6h with a* and the monotonic 
evolution of Lj^t with Ojj as shown in Figs. 2 and 4, respec- 
tively. 



2.2. Intrinsic time dependence of jet luminosity 

Based on the BH evolution we can derive the time profile 
of Ljet- Considering that the angular momentum is trans- 
ferred from the rapidly rotating BH to the disk, on which a 
positive torque is exerted, we think that the accretion onto 
the BH is probably halted. This state is discussed as the 
suspend accretion state by van Putten & Ostriker (2001), 
and as the nonaccretion solution by Li (2002). In this state, 
the BH evolution is governed by the BZ and MC processes. 
Based on the conservation of energy and angular momen- 
tum we have the evolution equations of the BH as follows. 



dM/dt = -{PBZ + PAIc), 



da,/dt = -M-^{Tbz +Tmc) 
= Bj^MA{a^,n), 



and 



BZ 



(9) 



Pmc) 
(10) 



A{a,,n) = -{Tbz +Tmc)/To + 2a,{PBz + PMc)/Po-{n) 

The powers and torques for the BZ and MC processes 
are expressed as (Wang et al. 2003), 



Pbz = Pbz/Po 



2ai 



fc (1 - k) sin^ 



2-{l-q) sin^ e 



Tbz = Tbz /To = 4a, (1 -I- q) 



(1 - fc)sin^ 
2-{l-q) sin^ 



MC 



MC I 



2ai 



'^/^ /3(1 -/3) sin^ 



2~{l-q) sin^ 



T. 



MC 



MC 



'To = 4a, (1 + 9) 



7r/2 



(1 - ^) sin^ 9d 
2 - (1 - g) sin^ 



(12) 



(13) 



(14) 



(15) 



where Tq = Bj^M^, and (3 = Hd/^h is the ratio of the 
angular velocity of the disk to that of the BH. 

Substituting Eqs. dH) and ^ into Eq. dH), we have the 
curves of Ljet{t) varying with time t for different values of 
'0 and n as shown in Fig. 6. The magnetic field on the BH 
horizon Br = lO^^G, the initial BH mass M (0) = TM© and 
the initial BH spin a,(0) — 1 are assumed in calculations. 

As shown in Fig. 6, the time profiles of Ljet exhibit a fea- 
ture of fast-rise and slow-decay, which arises from the non- 
monotonic evolution characteristic of 6h in terms of the BH 
spin as shown in Fig. 2. Thus the intrinsic time profile of the 
jet luminosity is derived naturally in our model rather than 
assumed ad hoc, and it resembles that given in PZLL99. 
The shapes of the profiles can be adjusted by changing the 
values of n and ip. It has been argued in LWM05 that the 
durations of various GRBs can be fitted by changing the 
value of n for the given initial BH spin. 
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Fig. 6. The curves of Lj^t vs time t for different n and 
observed angle "0, each curve is normahzed to unity at its 
maximum: a) for ■0 = 1° with n =5, 4.5, 4 and 3.5 in soUd, 
dash-dotted, dashed and dotted hnes, respectively; b) for 
n = 4 with V = 1°, 2°, 3° and 4° in solid, dash-dotted, 
dashed and dotted lines, respectively. 



3. Fitting observed light curves 

3. 1 . P recessing jet 

The misalignment in the spin axis of the BH and the angu- 
lar momentum axis of the binary causes the accretion disk 
around the BH to precess, which is known as the slaved 
disk precession. PZLL99 investigated the precession of the 
accretion disk in the gamma-ray binary, in which a neu- 
tron star fills its Roche lobe and transfers its mass to a 
BH. It is shown that this system has a precession period 
of about a second. The Lense-Thirring precession appears 
for a Kerr BH, if the accretion disk is inclined with respect 
to the equatorial plane of the BH (Lense & Thirring 1918). 
Recently, Reynoso et al. (2006, hereafter RRS06) studied 
the possible effect of Lense-Thirring precession on neutrino- 
cooled accretion disks. It is found that the precession period 
can be much less than 1 second. 

As the magnetic field is assumed to be anchored in the 
disk, this would result in the jet perpendicular to the mid- 



plane of the disk, and the precession and nutation of the 
disk directly implies the precession and nutation of the jet. 
To describe the jet kinematically we adopt the angular evo- 
lution of the jet in terms of the spherical angles 9 jet and 
(fijet given in PZLL99 as follows, 

4>jet = ^pre{t - to) + sin(il„„t), (16) 

Ojet = 0%t + §^ tan 0%, cos(l]„„i), (17) 

where ilpre = '^T^/Tpre ^ud flnu — ±27r/r„„ are the an- 
gular frequencies of precession and nutation with periods 
Tpre and Tnu, respectively. The direction of the observer 
can be assigned as Oobs and (pobs- A schematic picture of 
the precessing disk and jet is shown in Fig. 7, and the an- 
gle between the observer and the central locus of the jet 
is 

^ = COS~'^ {fobs ■ r jet) (18) 




Fig. 7. A schematic picture of a precessing disk with jet 
(adapted from RRS06). 



3.2. Fittings 

Substituting Eq. (HH) into Eq. we obtain the gamma- 
ray luminosity of a prcessing jet Ljet{t), which is adopted to 
fit the observed complex light curves of GRBs. The model 
contains six free parameters to fit the light curves (as de- 
scribed in Table 1), and the fitting procedures are described 
as follows. 

The first step is to determine the background. Following 
PZLL99, this is done by averaging the count rate of the 
initial, ~1800 data. The background is subtracted in our 
fitting process. 

Next, we introduce the start time Tstart and the end 
time Tend, which are the times of the start-up and shut-off 
of the real burst, respectively. As a rough estimate, we take 
Tstart and Tend as the first and last time to reach about 1% 
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Table 1. Free parameters that may vary per burst 



Parameters 


Note 


n 


The power-law index for the variation 




of the magnetic field on the disc 


'Tpre 


Precession period 




Frequency ratio of nutation to precession 


"jet 


Initial angular coordinate of a jet 


dobs, 4>obs 


Observation angles 



of the peak count rate. Our model fits the hght curve from 

^start to T^^ifi. 

From Figs. 2 and 6 (a), we find that the duration of 
GRB is sensitive to the value of n. Therefore, the first 
free parameter n is chosen to satisfy the observed dura- 
tion {Tend — Tgtart) in our model. Substituting n into Eqs. 
pUj) and pT|) , we obtain the evolution functions of BH mass 
M{t), BH spin a^{t) and half-opening angle 9 nit). The pre- 
cessing effect is described by the function ip{t). Substituting 
these functions into Eq. ([5]), we can produce the time vari- 
ability of the observed flux. 

Finally, we use the simulated annealing to fit the ob- 
served bursts. The other five parameters in Table 1, i.e.,rpre, 
Rn, Q'jet^ ^obs and (f>obs are chosen freely in this simulation. 
After each iteration we determine the from the fitting 
of Ljet{t) to the observed burst profile. This value is mini- 
mized by the annealing algorithm. 

We apply this model to fit several observed GRBs by ad- 
justing the value of the above six free parameters (see Table 
1). The data of the light curves are taken from BATSE and 
HETE. We fit the light curves of GRBs from the third en- 
ergy channel, i.e., 100-300keV for BATSE, and 30-400keV 
for HETE. The fittings results are presented in Fig. 8 with 
the parameters listed in Table 2. 

Inspecting Fig. 8, we find that the light curves of these 
GRBs can be fitted by our model. Some simple light curves 
of GRBs are fitted in a rather satisfactory way as shown 
in Figs. 8a, 8c, 8e and 8g, and complicated light curves are 
fitted with good overall shapes as shown in Figs. 8b and 
8d. Fig. 8f shows the fitting for the more complicated burst 
GRB 030519a for which most of the peaks are well fitted 
except for the third one. GRB 920701 is a peculiar asym- 
metric burst (Romero et al. 1999; RRS06), with a slower 
rise and a faster decay compared with the other bursts. It 
turns out that this type of burst can be also fitted by our 
model as shown in Fig. 8b. 

4. Discussion 

In this paper we discuss a model for fitting the light curves 
of GRBs based on the coexistence of the BZ and MC pro- 
cesses, combining the evolution of an extreme Kerr BH sur- 
rounded by a precessing disk. Comparing with PZLL99 , 
PZTOl and RRS06, our model can explain both the intrin- 
sic time profile and the observed duration of GRB. The 
fittings are done with only six parameters. However, there 
are still several issues related to this model and the fittings. 

First, a precessing jet tends to produce light curves 
in which subsequent peaks appear at periodic intervals. 
However, the observed gamma-ray bursts show no evidence 
for periodicities at any time-scale. PZTOl solved this prob- 
lem by considering the interaction between the jet and the 
interstellar medium. In our model, the distance from the 



BH to the emission region is fixed. If this distance varies 
during the GRB, the asymmetry of the subsequent peaks 
in the light curves could be explained by our model. 

Second, the individual peaks produced by precession 
clearly lack the observed strong asymmetries, i.e., a fast rise 
and a slow decay (Fenimore, Madras & Nayakshin 1996). 
In our model, the time variability of the central engine is 
rather smooth. Thus our model cannot explain the strong 
asymmetries. As described in PZTOl, the effect of the cur- 
vature of the jet front or the cooling of gamma-ray emitting 
electrons can make the pulse profile a 'fast rise and slow de- 
cay'. We will take these effects into account in future work. 

Third, this model is good in fitting the overall shapes 
of the light curves of the bursts, but it is difficult to fit 
the smaller details such as the variability on time scales 
of milliseconds. Recently, Wang et al. (2006) discussed the 
modulations of the screw instability in the BZ process on 
the light curves of GRBs, and they argued that the variabil- 
ity timescales of ten milliseconds can be interpreted by two 
successive flares due to the screw instability of the magnetic 
field. It is also shown that the individual peaks produced 
by releasing and recovering magnetic energy have a fast rise 
and slow decay profile. We shall improve the fitting of the 
light curves of GRBs by combining the jet precession with 
screw instability of the magnetic field in future work. 

We cannot perform satisfactory fits for all kinds of GRB 
light curves with this simple model, especially for those dis- 
playing extraordinary periods of quiescence between bursts. 
Fig. 8f is one of these cases, where a time-lag in the third 
peak is not reproduced. This type of burst is hard to fit 
by the jet-time-profile as given in Fig. 6, i.e., the feature of 
fast-rise and slow-decay. The combinations of parameters 
listed in Table 2 are not unique in the fittings. For exam- 
ple, the light curve of GRB 920701 can be fitted as well as 
Fig. 8b with the parameters (n=3.10; Tpre— 5.18, i?o=8.36; 
0%t^l.O2, 0ofc,=1.19, 0ofc,=242.58). ' 
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Fig. 8. Fits of the light curves from the third energy channel of (a) GRB 910717, (b) GRB 920701, (c) GRB 990123, (d) 
GRB 001225, (e) GRB 030329, (f) GRB 030519a and (g) GRB 031111a. The real bursts profile and the fitting curves 
are plotted in full and broken curves, respectively. 

Table 2. Values of the parameters for fitting several GRBs 



GRB 


Tstart (s) 


Tend (s) 


n 


'Tpre (s) 


Rn 


e'^et (deg) 


dobs (deg) 


(j)obs (deg) 


x' 


910717 


1.86 


7.94 


3.15 


7.00 


4.22 


1.78 


0.89 


39.16 


0.88 


920701 


1.02 


6.98 


3.10 


7.24 


20.01 


2.63 


0.69 


188.25 


1.52 


990123 


0.51 


100.86 


4.10 


128.37 


1.11 


3.89 


2.97 


149.17 


2.05 


001225 


0.28 


33.37 


3.60 


14.68 


6.08 


0.74 


4.28 


172.58 


0.84 


030329 


7.24 


51.99 


3.70 


70.89 


2.60 


11.16 


11.20 


69.88 


1.78 


030519a 


0.09 


14.18 


3.26 


7.80 


4.35 


1.54 


0.41 


259.86 


1.63 


031111a 


2.56 


16.16 


3.26 


22.13 


12.86 


5.58 


4.96 


101.63 


0.14 



